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O . ABSTRACT 
CN 

O I Context. The spectra of massive binaries may be affected by interactions between the stars in the system. These are believed to pro- 

^ , duce observational phenomena such as the Struve-Sahade effect. 

■ Aims. We simulate the spectra of massive binaries at different phases of the orbital cycle, accounting for the gravitational influence of 

the companion star on the shape and physical properties of the stellar surface. 

Methods. We used the Roche potential modified to account for radiation pressure to compute the stellar surface of close circular 
systems. We further more used the tidal interactions with dissipation of energy through shear code for surface computation of eccen- 
tric systems. In both cases, we accounted for gravity darkening and mutual heating generated by irradiation to compute the surface 
temperature. We then interpolated non-local thermodynamic equilibrium (NLTE) plane-parallel atmosphere model spectra in a grid to 
obtain the local spectrum at each surface point. We finally summed all contributions, accounting for the Doppler shift, limb-darkening, 
and visibility to obtain the total synthetic spectrum. We computed different orbital phases and different sets of physical and orbital 
(— I . parameters. 

Results. Our models predict line strength variations through the orbital cycle, but fail to completely reproduce the Struve-Sahade 
I I effect. Including radiation pressure allows us to reproduce a surface temperature distribution that is consistent with observations of 

O ■ semi-detached binary systems. 

Conclusions. Radiation pressure effects on the stellar surface are weak in (over)contact binaries and well-detached systems but 
C/J ' can become very significant in semi-detached systems. The classical von Zeipel theorem is sufficient for the spectral computation. 

' Broad-band light curves derived from the spectral computation are different from those computed with a model in which the stel- 

lar surfaces are equipotentials of the Roche potential scaled by the instantaneous orbital separation. In many cases, the fit of two 
Gaussian/Lorentzian profiles fails to properly measure the equivalent width of the lines and leads to apparent variations that could 
explain some of the effects reported in the literature. 

Key words. Stars: massive - Binaries: general - Stars: fundamental parameters - Stars: atmospheres - Binaries: spectroscopic 

O 

(N 

. 1 . Introduction The first models of binary systems in which the stars are not 
. spherical have been proposed by Russell & Merrill (11952) to 
O ■ Studing massive stai's is important because they play a key role reproduce the light curves of binaries. Then, Kopal (fT959] l in- 
rn ; m galaxy evolution: their strong winds interact with the ambi- Produced the Roche potential approach that was first used by 
ent interstellar medium and can trigger the formation of new Lucy (1968 ) and Wilson & Devinney (TWT) followed by nu- 
stars, they produce heavy chemical elements, and they are an ^^^.^^^^ ^ther works. In a first step, we designed a similar code 
important source of UV radiation. Recent research indicates that ^j^^^^j^. ^^^^^^^ ^inaiy systems (Palate & Rauw|20l2] here- 
a large part of massivestai-s ( ~ 50%) form abinai-y or a multiple ^^g^ p^p^^ j) j„ ^hich the distorted shapes of the stellar surfaces 
H ■ system (Mahy et al. |2009l Sana & Evans|20lD- The specti-a of ^^^^g^ ^y gravitational interaction. The emerging spectra are 
y . binary systems are an invaluable source of information for de- ^^^^ computed accounting for the shape of the stars and their lo- 
termining the physical properties (such as masses, temperatures, properties. In this paper we extend the method to incoiporate 
and radii) of stars. It has been shown, however, that the obser- radiation pressure effects on the stellar surface and the case of ec- 
vational analysis of massive bmanes is complicated by effects ^^^^^^ Uneivies. This new version of our CoMBiSpeC (code of 
that are linked to interactions in these systems tha^ffect the ^^^^^^^ binary spectral computation) model now allows spectral 
spectra and the specti-al classification (Sana et al. 120051 Linder computation of (almost) any massive binary system, 
et al. 2007 Linder 2008 and references therein). Unfortunately, 
the majority of spectral modelling codes are designed for iso- 
lated spherical stars and cannot reproduce the particular effects 2 Improvements Of the method 
of the binarity. Therefore, it is important to improve our mod- 
els of spectral computation to accurately represent the spectra Our method for modelling the spectra of circular binary systems, 
of these stars. In this context, we have developed a method that presented in Paper I, consists of computing the actual stellar 
takes into account some of the effects that are produced by grav- surface following the Roche potential approach and accounting 
itational interactions in the system. for the local gravity and temperature on each surface element. 
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Non-L TE spectral model grids (TLUSTY, Lanz & Hubenv l2003l 
12007 ') are used to compute the integrated spectrum of the star at 
each orbital phase (see also, e.g., Linnell & Hubeny [19941 and 
Linnell et al. 2012). Because our spectra are synthetic, they are 
free of observational noise. In this section we describe the two 
improvements that have been implemented to the method. The 
first is the inclusion of radiation pressure effects; i.e., the irra- 
diation of each star by the companion and the effect of the ra- 
diation on each star's own surface. The second is the use of an 
explicit calculation of the stellar surface, from first principles, 
which here is applied to the case of eccentric binaries. In ec- 
centric systems, the orbital separation changes as a function of 
phase, and accordingly, the deformations of the surface and the 
orientation of the tidal bulge are also variable. 

2.1. Radiation pressure effects 

There exists quite an extensive body of literature (see e.g. 
Pen nine et al. |2009l Howarth [T997al Sch uerma n [T972] Drechsel 
et al. [T995] and Phillips & Podsiadlowski'200?) on the impact of 
radiation pressure on the Roche potential in early-type binaries. 
According to Dermine et al. (12009) . the effects of radiation pres- 
sure in a binary system can be separated into three parts: radia- 
tion pressure on each star's own surface, the effect on the com- 
panion, and radiation pressure on the matter outside the binary 
orbit. Here we focus on the two first effects, which are called in- 
ternal and external radiation pressure (Drechsel et al. I1995I ). The 
impact of the radiation pressure on the matter outside the binary 
is, of course, very important for the surrounding medium but is 
irrelevant for calculating the shape of the stars. Howarth ( l997al ) 
pointed out that the effect of the internal radiation pressure can 
be treated as a simple scaling of the Roche potential of the stars. 
Another approach, proposed first by Schuerman C1972\ consists 
of scaling the mass of the stars rather than the entire potential. 
The former approach is more appropriate because the emitted 
radiation bolometric flux and hence the radiation pressure scale 
with the local gradient of the potential (which also includes the 
attraction by the other star and the centrifugal forces). Therefore, 
according to Howarth (1997a), the stellar topology is not affected 
by internal radiation pressure and the latter simply scales the lo- 
cal surfaces gravity. The Newtonian gravity is defined by 



g 



-V0, 



where is a scalar gravity potential per unit mass. 
The von Zeipel (119241) theorem can be writteifl 



(1) 



(2) 



The magnitude of the radiative acceleration is given by 
(Howarfli d997al ) 

Cmd = T, (3) 

c 

where k is the flux mean opacity per unit mass and T' is the 
astrophysical bolometric flux (f oc T'^). 

If we define F as the ratio of the radiative to the gravitational 
acceleration, Oiad = ^g- According to Howarth (1997a), this ratio 
is constant for a given star and independent of the position on the 
stellar surface since both accelerations vary with T"*. Therefore, 
we have that the effective potential is the potential without radi- 
ation pressure minus the radiative acceleration 



geff = (1 -r)g. 



(4) 



' Throughout this paper, the temperature is to be understood as effec- 
tive temperature. 



From equation [T] (pgff = (1 - r)(p and the ratio F can be 
written as (following the Castor et al. ll975l theorv for evaluating 



r = -rig 

c 

O" Th ^4 

= crT 

MhC 



(5) 



with ^ w 0.036 m^kg crj-/, the Thomson-scattering cross 
section and cr the Stefan-Boltzmann constant. 

Because F is independent of the position on the stellar sur- 
face, we can compute it for T - Tpoie and g - gpoie 



cr Th 



1 



\Spole\ 



(6) 



The advantage of this method is that the von Zeipel theorem 
remains valid. Another advantage is the simplicity of the method 
which does not require many computational resources. 

The external radiation pressure is more difficult to treat. We 
have based our method on the approach of Drechsel et al. (119951 ) 
and Phillips & Podsiadlowski (2002). This method consists of 
scaling the mass of the companion in the Roche potential. The 
scale parameter 6 - is computed iteratively for each surface 

point. The external radiation pressure can be seen as a force that 
decreases the attraction of the companion. Therefore, the com- 
panion seems less massive, and consequently the mass has to be 
scaled. The new "Roche" potential can be written 



where r 



q(l-6(r,0,<p)) 



q + I 



2rcos(^sin0 + 1 



■ sin^ 
2 (7) 

-qr cos (fi sin , 

- rcos(fsm0, y = 
rsini^sinfl, and z = rcos0. Here, and (f are the colatitude 
and longitude angle in the spherical coordinates centred on the 
star = and ^ = correspond respectively to the north pole 
of the star and the direction towards its companion. 
The radiation pressure can be written 
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-\- -\- z , q 



rad ' 
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(8) 



where "Pi is the angle between a local surface normal on the 
irradiated star and the direction to a surface element on the other 
star ly is the specific intensity of the external radiation field, and 
V is the frequency. 

The radiation pressure gradient per unit mass is equal to (in 
plane-parallel approximation) 



rad 



_ldP,-ad 
P 



dr 



dly 



= 7l/o°°'^v^cos2»I'idydw 
= ^I/o"'^v/v cos^-idydw. 



(9) 



Following Castor et al. (119751 , we assume that the dominant 
source of opacity comes from electron scattering, thus, we have 

Ky = ^. 

The value of 6 is therefore computed with an iterative 
method. We start from the classical unmodified Roche potential 
and compute the radiation received by each surface element from 
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its companion. The radiation is computed by summing the con- 
tribution of all points that are visible from the primary (resp. sec- 
ondary) on the secondary (resp. primary). The method applied 
to compute the received bolometric flux by a point is derived 
from equation |9] and is similar to the reflection effect treatment 
of Wilson dl 9901 1 ■ Indeed, we have ly - ly^ x (I - u + u cos ^'2), 
where ^'2 is the angle between a local surface normal on the 
emitting star and the direction towards an irradiated surface el- 
ement on the other star, and u is the linear limb-darkening co- 
efficient based on the tabulation of Claret & Bloemen (1201 11 1. 
Moreover, we can write doj - ^^2££^h. where s is the distance 

s- 

between the irradiated and emitting surface. Finally, if the radia- 
tion comes only from the companion, we have lyfi dv = 

-, where d is the 



CM 



The local 6 factor is thus defined hy S - 
distance between the point and the centre of the emitting star. 
This factor can be written (for the radiation emitted by the sec- 
ondary and received by the primary) as 



crcr,i,d 
ncmHGM2 



^ r^(l -M + MC0S^2)C0S>PiC0S*P2^ 



co.s4'i>0 
co.s4'2>0 



(10) 

Then, with the S factor, we recompute the stellar surface with 
the modified Roche potential. Again, we evaluate the 6 factor 
and surface until we reach convergence. 

2.2. Eccentricity 

The extension to eccentric systems uses the TIDES0 code 
(Moreno et al. [T999ll2005ll20TTT l. This code computes the time- 
dependant shape of the stellar surface for eccentric and/or asyn- 
chronous systems. It also provides velocity corrections for rigid 
body rotation. The method consists of modelling a deformable 
and perturbed surface layer that lies upon a rigid body that is 
in uniform rotation by solving the equations of motion of a grid 
of surface elements. These equations take into account centrifu- 
gal and Coriolis forces, gas pressure, and viscous effects in ad- 
dition to those of the gravitational potential of both stars. It is 
important to note that in asynchronously rotating or eccentric 
systems, the linear approximation of the Roche potential is no 
longer valid and the presence of viscosity in the medium may 
lead to non-linear effect^ The surface shape and local veloci- 
ties are obtained as a function of orbital phase from the solution 
of the equations of motion. Full details of the model are given 
in Moreno & Koenigsberger (1 19991 1. Toledano et al. (2007), and 
Moreno et al. (2005,2011). 

With the derived values of radius and velocity perturbations, 
we compute the surface gravity and temperature. Assuming the 
perturbations to be small, the gravity is computed by using a 
classical Roche gradient. 

The temperature distribution is computed following the von 
Zeipel theorem 



l|Vi2iocal 



local 



- pole 



VQ, 



pole ^ 



(11) 



where the gravity-darkening parameter p - I'm the case of mas- 
sive stars. 

We also include reflection effects that heat the surface of 
the star through irradiation of each component by the other. 



^ Tidal interactions with dissipation of energy through shear. 
^ The viscous stresses are included in the equations of motion and 
connect the surface elements to each other and to the inner rigid body. 



The effect was treated by following the approach of Wilson 
(119901) . We neglected any cross-talk between the surface ele- 
ments (e.g. due to horizontal advection or radiative exchanges 
between neighbouring surface elements of different tempera- 
ture). This is probably an important approximation because this 
effect might smooth out some of the temperature variations, es- 
pecially near periastron (for eccentric systems) where the stars 
move faster and where the gravitational interactions are stronger 
These interactions and velocity imply that energy dissipation due 
to viscosity effects is more significant near periastron and thus 
the cross-talk is expected to be higher 

The spectral computation was performed with the 
CoMBiSpeC code, which was described in Paper I. For 
this paper, some minor changes have been applied to take into 
account particularities linked to the eccentricity, such as the 
variation of the separation between the two stars. 



3. Circular-orbit models 

We present here models of four binary systems with circular 
orbits. Three of them were previously studied in Paper I, i.e., 
models 1, 2, and 3. These models are based on the HD 159176, 
HP 16 5052, and HD 100213 systems analysed by Linder et al. 
( 120071 ). The fourth system, model 4 in this paper, is based on Sk- 
67° 105 studied by Bonanos (2009 ). All four of these systems are 
in synchronous rotation and, thus, are in equilibrium configura- 
tion so that the Roche potential approach is valid. For all mod- 
els, we compared the shape, temperature, and gravity variations 
caused by radiation pressure. Table[T]gives a summary of the pa- 
rameters (from Linder et al. 12007 1 and Bonanos 12009) used for 
the computation. We have chosen the polar temperature (rather 
than a surface mean temperature) as the reference effective tem- 
perature in Table [T] This is because Tpoie is used to compute the 
local temperature via the von Zeipel theorem and remains essen- 
tially unaffected by radiation pressure effects. Finally, for model 

4, we investigated the Struve-Sahade effect for some lines. 

3.1. Models 1 and 2 

We have analysed the spectr^ of models 1 and 2 at 20 orbital 
phases and found that the impact of radiation pressure is rather 
low. The comparison was made by subtracting the normalized 
spectra with radiation pressure to the spectra without radiation 
pressure. On each of these "residual" spectra, we searched for 
the maximum, mean, and median values and the standard devi- 
ation. Then we searched for the maximum of the maxima (over 
the orbital phase) and computed a mean of the medians, the stan- 
dard deviations, and means. We also looked for the spectral bins 
that present deviations higher than 1 % (over a total of 360000 
bins). Table|2]gives a summary of these values for the combined 
spectra of the studied models. For example, the highest differ- 
ences are 0.01 13 and 0.0191 for models 1 and 2. These maxima 
occur in both cases at phase 0.5, i.e. at the conjunction with the 
primary in front of the secondary. These differences are signif- 
icant but, the mear0 difference is only ~7 - 9 x 10 which is 
very small. The general appearance of the spectra is preserved 
and the conclusions reached in Paper I remain valid. We can see 
in Fig.[T]that the spectra of model 1 with and without radiation 
pressure are nearly indistinguishable (top panels). 



Throughout this paper, the word "spectrum" refers to a synthetic 
spectrum. 

' We only took into account the spectral bins that include spectral 
lines for this mean value. 
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Table 1. Parameters (from Linder et al. 120071 models 1 to 3, and Bonanos l2009l model 4) of the circular binary systems simulated 
in this paper Inclinations in brackets stand for non-echpsing systems. 



Parameters 


Model 1 


Model 2 


Model 3 


Model 4 


Period (day) 


3.36673 


2.95515 


1.3872 


4.251 


Mass ratio 


0.96 


0.87 


0.68 


0.42 


Semi-major axis (Rq) 


38.23 


31.25 


17.34 


38.94 


Inclination (°) 


(48) 


(23) 


77.8 


89.9 


Mass of primary (Mq) 


33.8 


25.15 


21.7 


30.9 


Mass of secondary (Mq) 


32.41 


21.79 


14.7 


13.0 


Primary polar temperature (K) 


38000 


35500 


35100 


35000 


Secondary polar temperature (K) 


38000 


34400 


31500 


32500 


Polar radius of primary (Rq) 


9.37 


9.10 


6.74 


15.1 


Polar radius of secondary (Rq) 


8.94 


8.47 


5.62 


11.14 



The largest and mean radii of the stars are smaller if radi- 
ation pressure effects are included. The radiation pressure also 
decreases the mean \og(g) by ~ 0.05 despite the smaller radii. 
Finally, the mean temperatures increase by lOOK. This increase 
of the temperature and decrease of gravity cause a weak increase 
of a few per cent of the total flux emitted by the stars. Fig. [T] 
lower panel, displays the log of the ratio of the fluxed spectra 
with and without radiation pressure and underlines the small dif- 
ferences in the continuum. Table [3]gives a summary of the radii 
and temperatures of circular models. 



Table 3. Temperature and radius of stars at different locations at 
the stellar surface. 



Stars 




^point 


Rbxi. 


^side 




T'back 










(K) 


(K) 


(K) 


Model 1' 


9.51 


9.46 


9.70 


37442 


38061 


36705 


Model 1- 


9.06 


8.97 


9.23 


37496 


38271 


36800 


Model 2' 


9.32 


9.44 


9.62 


34654 


34937 


33577 


Model 2- 


8.66 


8.67 


8.96 


33643 


34465 


32510 


Model 3' 


7.12 


8.44 


7.63 


33198 


29306 


29211 


Model 3' 


5.87 


6.92 


6.43 


30112 


30855 


27369 


Model 4' 


15.80 


15.95 


16.42 


33411 


33882 


32121 


Model 4- 


11.61 


10.72 


12.89 


31185 


34062 


27912 



Notes. primary star, -: secondary star, side: r(^,6) = r(n/2,n/2), 
point: r((p,e) = r(0,n/2), back: r(4),0) = r(n,n/2) 



3.2. Model 3 

This model leads to an (over)contact binary if radiation pressure 
is not included, but when it is included, the system becomes de- 
tached0. Surprisingly, however, the impact of radiation pressure 
on the spectrum is weak. This is because the parts of the star 
that contribute most to the spectrum come from the rear and side 
parts of the star which are less modified by the radiation pressure 
because S is small (see Fig|2]l. Moreover, the shape of the star is 
not strongly affected by the radiation of the companion. This is a 
priori unexpected but the reason for it is that the most deformed 
parts of the stars mainly see the coldest part of the companion, 
which leads to a rather low value of the maximum S. We stress 
here that in systems where the two stars are in (over)contact (or 
nearly so) and are of similar size and temperature, the reflection 



* Strictly speaking, the CoMBiSpeC model cannot handle overcon- 
tact configuration, but the most extreme cases correspond to situations 
where both stars fill their Roche lobe. 




4000 4100 4200 4300 4400 4500 4600 

Wavelength (A) 

Fig. 1. Top: Comparison between the normalized spectra com- 
puted without radiation pressure (in black) and with radiation 
pressure (in grey) for model 1. To better distinguish spectra, 
the black line is wider than the grey line. Orbital phases are, 
from bottom to top: phases 0.1, 0.25 and 0.4. The spectra are 
shifted vertically by 0.25 continuum units for clarity. Middle: 
Corresponding residuals (spectrum with radiation pressure mi- 
nus spectrum without radiation pressure). The residual plots are 
shifted by 0.015 for clarity. Bottom: Logarithm of the ratio of the 
fluxed spectra with radiation pressure divided by the one with- 
out. These plot are shifted by 0.006 for clarity. 



effect is not sufficient to counterbalance the gravity darkening. 
Thus, the coldest parts of the stars are those near the LI point. 

In Paper I, we pointed out that we failed to reproduce the 
surface temperature distribution deduced by Linder et al. (2007)) 
for HD 100213. These authors showed in thek analysis through 
radial velocity measurements that the He II and He I lines were 
not formed in the same region of the stellar surfaces. The small 
changes of the stellar surface induced by the radiation pressure 
increase the importance of the reflection process and are suf- 
ficient to now explain the surface temperature distribution ob- 
served by Linder et al. The He II lines have a smaller radial 
velocity amplitude and thus are apparently stronger in the hemi- 
sphere facing the companion, whilst the He I lines have larger 
velocity amplitudes, indicating that they are stronger on the op- 
posite hemisphere. This suggests that the hemisphere facing the 



4 



M. Palate et al.: Spectral modelling of massive binary systems 
Table 2. Comparison between the spectra with and without radiation pressure. 



Models 


l|A„,„,|| 




ll^/)ie<Y/(3« II 


IIo-aII 


A > 1% 


Model 1 


0.0113 


6.8 X 10-^ 


3.2 X 10-* 


0.0014 


12 


Model 2 


0.0191 


8.7 X lO^* 


3.8 X 10-* 


0.0018 


213 


Model 3 


0.0063 


5.1 X lO-'* 


3.0 X 10-* 


8.7 X lO-* 





Model 4 


0.1247 


0.0069 


0.0022 


0.0171 


95820 



companion is hotter than the opposite one. When radiation pres- 
sure is added, the stars are less deformed. This implies a weak 
temperature increase of the facing hemispheres caused by the 
decreasing gravity darkening and a higher effectiveness of the 
reflection process. This weak increase is therefore sufficient to 
explain the observations. 




4 8 12 

X (fle) 



16 20 



I 



0.03 

0.025 

0.02 

0.015 

0.01 

0.005 





Fig. 2. Value of the radiation pressure parameter, 5, over the stel- 
lar surface for model 3. The highest values are located near the 
LI point. 




0, 



-20, 
20 

y(^0). 




Espinosa Lara & Rieutor d (120 IH and Maeder (11999) have 
suggested that von Zeipel's (119241 ) classical theorem requires 
modifications. These theoretical studies present more complex 
models that could be, at first approximation, equivalent to the 
classical von Zeipel theorem with a gravity-darkening parame- 
ter (GDP) p smaller than 1 .00. Therefore, we chose this very de- 
formed system in which the von Zeipel theorem implies strong 
temperature vaiiations to test the influence of p on the spectra. 
We computed the spectrum for values of the GDP of 0.25, 0.50, 
0.75, 1.25 and 1.50, all other parameters having the same value 
as model 3 in Table [1] We found that for a moderate change of 
p (of +/ - 0.25) the spectra are not notably affected. We also 
found that the variations remain within the noise level of real 
data for variations of up to +/ - 0.50. Therefore, since we made 
some other impoitant assumptions (no cross-talk, assumptions 
underlying the Roche potential), it is not useful to adopt another 
more complex formalism for gravity darkening because it does 
not notably affect our simulations. 

3.3. Model 4 

This system was chosen as an example of one in which the radi- 
ation pressure considerably modifies the shape of the stars (see 
Fig. |3]l, underlining the potential significance of this effect. The 
system parameters (see Table [1} are based on the observational 
study of Sk-67°105 that was carried out by Bonanos (.2009 « . 



X (Rq) 

Fig. 3. Surface modification by the radiation pressure for model 
4. Top: Star surface without radiation pressure effect. Bottom: 
Star surface with radiation pressure effect. 

3.3.1. Spectral classification 

Observational studies have shown that the primary and sec- 
ondary are an 07-8V and an 08-8.5 III-V star. We first deter- 
mined the spectral type of our synthetic spectra with the Conti 
(Conti & Alschuler 1971) - Mathys ( 1988 1989) criterion (here- 
after Conti-Mathys ciiterion). This criteiion is based on the ratio 
of the equivalent widths (EWs) of He I ^4471 to Hell ^4542 
for the spectral-type determination and on the ratio of the EWs 
of Si IV A 4089 to He I /1 4143 for the luminosity-class determi- 
nation. We found that our model spectra that include radiation 
pressure correspond to an 06.5-7I type for the primary spectrum 
and an 07. 5-081 type for the secondary. The determination of 
the spectral type gives the same result whether the measurement 
is made on individual spectra or on combined spectreQ. In our 
model, the stars seem hotter, and the problem encountered in 
Paper I concerning the luminosity class is again present because 

' "Combined spectrum" refers to the spectrum of the primary plus 
the spectrum of the secondary, i.e. the spectrum of the entire system. 
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we overestimate the luminosity. In a second approach, we used 
the Walborn & Fitzpatrick (1990^ atlas to compare individual 
spectra to the reference atlas and thus re-classify the stars. With 
the atlas we found for the primary and secondary an 0(6.5)-7 V 
and an 08-8.5 (III)-V classification. The classification using the 
atlas is of course more qualitative, therefore we give a range and 
the best correspondence is in brackets. Within the uncertainties, 
the agreement with the observations is good. 

3.3.2. Struve-Sahade effect 

The Struve-Sahade effect (hereafter S-S effect) was originally 
defined as the apparent strengthening of the secondary spectrum 
when the star is approaching the observer and its weakening as it 
moves away (Howarth et al. 1997b, Linder et al. 2007 1. A more 
general definition is the apparent variation of the line strengths 
of either of the binary components as a function of the orbital 
phase (Howarth et al. 1997bi) . Bonanos ('2009') reported a quite 
considerable S-S effect for several spectral lines of Sk-67°105. 
We studied eight of these Hnes in detail: He I AA 4026, 4143, 
447 1 , 47 1 3, 50 1 6, Si IV ^ 4089, He II 4542, and 54 1 1 . 

First, we analysed the EWs on individual spectra of the pri- 
mary and the secondarjQ. The measurements of the line EWs 
were performed with the MIDAS software developed by ESO. 
The EWs of lines in the spectra of the individual stars were de- 
termined directly by simple integration with the integrate/line 
MIDAS routine, whilst for the binary system, we used the de- 
blend/line command as we would do for actual observations of 
a real binary spectrum. The latter routine fits two Gaussian line 
profiles to the blend of the primary and secondary lines. 

The EWs of the different lines measured on the individ- 
ual spectra display phase-locked variations. For the primary, the 
EWs of the He I lines reach a maximum at phase 0.5 (when the 
primary eclipses the secondary) because the coldest part of the 
star is the rear part. Two minima are visible at phases 0.05 - 0.1 
and 0.9 - 0.95 because the observer starts to see the front part 
of the primary star. This part of the star is the hottest because 
of the very effective reflection in this system. The minimum is 
not observed at phase because of the roughly annular eclipse 
that hides a substantial part of the stellar surface. The variation 
amplitude is line-dependent with ~ 5% for He I /t 4026, 60% for 
He I /14143, and 20% for the other He I lines. As expected for 
heating effects, the He II lines exhibit a reverse variation com- 
pared to the He I lines (see FigHJl. The Si IV line displays a sim- 
ilar variation as the He I lines. 

For the secondary, the EWs display the same variation as 
the primary, but shifted by 0.5 in phase. We can notice the zero 
value at phase 0.5 when the secondary is totally eclipsed. The 
variation amplitude of the Hel ^^4026, 4713, 5016, and SilV 
lines is small and less than ~ 10%. The amplitude is greater than 
~ 20% for the other lines (we did not account for the zero value) . 

The measurements of the combined spectra using a two- 
Gaussian fit to the blended lines agree relatively well with the 
variations measured on individual spectra for several lines: Hel 
AA4026, 4143, 4713, 5016, Hell AA4542, and 5411. The vari- 
ation amplitude is sometimes overestimated, though. The EWs 
of Hel A 441 1 are asymmetric before and after phase 0.5. The 
variation does not seem regular but is similar to the variations 
measured on individual spectra, i.e. a decrease of the EW of the 



secondary and an increase for the primary during the first half of 
the orbital cycle. Finally, the Si IV line displays irregular varia- 
tions. 

Bonanos (2009) reported a strong S-S effect in Sk-67°105 
for several lines that we studied in the model. The variations ob- 
served are consistent with the definition of Linder et al. (I2007I I. 
i.e. in the first half of the orbital cycle the secondary (resp. pri- 
mary) line is deeper than the primary (resp. secondary) line and 
the situation is reversed in the second half of the orbital cycle. 
However, even though some lines in our synthetic spectra dis- 
play strong EW variations, we did not observe this type of phase 
dependence. 
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^ The mean relative error of the normalized EWs can be estimated to 
cr = 10"^ - 10"^. The cr given by the deblend/line routine is of the same 
order of magnitude. This error estimation is valid for all studied lines 
throughout this paper. 



Fig. 4. Example of EW variations of the Hell /1 4542 line. The 
EWs have been normalized to the highest value. Black: Primary. 
Grey: Secondary. Left: Variations measured on individual spec- 
tra (EWp, max = 0.53A, EWs, max = 0.47 A). Right: Variations 
measured on combined spectra by fitting two Gaussians to the 
blended hnes (EWp, max = 0.52A, EWs, max = 0.16A). 



3.3.3. Radial velocity and synthetic broad-band light curve 

The last two characteristics that we studied are the light curve 
in the wavelength range 3800 - 7100 A and the radial veloc- 
ity curve. The two eclipses are clearly visible in the light curve. 
The depth of the eclipses agrees well with the observations. We 
stress that light curves are "by-products" of the spectral compu- 
tation and that CoMBiSpeC is not primarily designed for light 
curve computation. The semi-amplitude of the radial velocity 
also agrees well with the observations. The radial velocity was 
computed by taking the mean value of the velocity of the visi- 
ble points of the stellar surface at a given phase weighted by the 
surface projected along the line of sight. The values are given 
for the primary and the secondary (observational value in brack- 
ets): 138.3 (137) km s"' and 332.9 (326) km s"' . Finally, the syn- 
thetic light curve and radial velocity curve agree well with those 
derived by Bonanos (I2009I I. 



4. Eccentric models 

We present here five models of eccentric binary systems. 
These models are inspired by the following massive bi- 
naries: HD93205, HD93403, HD 101131, HD 152218, and 
HD 152248, and are named models El, E2, E3, E4, and E5. 
The stars were studied by: Antokhina et al. (i2000 ). Rauw et 
al. (2000}, Gies et al. (|2UU2] |. Sana et al. (l2008l l. Mayer et al. 
(2008), and Sana et al. (2001 ). The parameters used for the com- 
putation are given in Table|4] We studied the spectral classifica- 
tion, the radial velocity curves, the broad-band light curve, and 
finally the S-S effect. We computed the surface of these stars 
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with the TIDES code (Tahle^bottom gives the specific param- 
eters used in the TIDES code) and used a modified version of 
the CoMBiSpeC code to compute the gravity and temperature. 
These quantities are not uniformly distributed across the stellar 
surface and, in addition, change over the orbital cycle. 

A sample of the phase-dependent behaviour of the mean ra- 
dius, gravity, and temperature as well as their maximum and 
minimum values are shown in Figures |5] |6] and Q The max- 
imum radius is largest near periastron (just before and/or just 
after). A model based on an instantaneous Roche lobe, such as 
the Nightfall^ pmgmm, presents the strongest variation strictly 
at the periastron passage. However, in the TIDES code, viscous 
effects that lead to a delay are taken into account. The most com- 
pact systems, model El, the secondary of model E2, and model 
E3, remain nearly spherical during the entire orbital cycle (see 
FigE) because the mean radii are nearly constant and the highest 
relative difference between the maximum and minimum radii is 
less than 0.5%. This is not surprising because models E2, and E3 
have the longest periods, hence the widest separation, and model 
El contains relatively compact main-sequence stars. 




Orbital phases 



Orbital phases 
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Fig. 5. Variations of the mean, highest and lowest visible ra- 
dius {left: the primary radii, right: the secondary radii). Top: 
Variations for model El. Bottom: Variations for model E3. 



Fig. 6. Variations of the mean, highest and lowest visible radius, 
temperature, and log(g) for model E4. Left: Values for the pri- 
mary star. Right: Values for the secondary star. 
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4.1. Spectral classification 

First, we appUed the quantitative Conti-Mathys criterion (Conti 
& Alschuler [T97n Mathvs [T988l[T989l l. The spectral types that 
we found for our simulated spectra either agree well with the 
observations or are hotter than the observations. The hotter stars 
present a weak He I /1 4143. This makes measuring the spectra 
of this system difficult and uncertain. The luminosity classes 
are, as in Paper I, overestimated for models El, E3, and E5. 
Model E2 and the primary of model El are too hot for the Conti- 
Mathys luminosity criterion. The luminosity classes agree well 
for model E4. The problem of simulating singlet transition of 
He I reported by Najarro et al. (2006 1 and discussed in Paper I 
is probably responsible for the general luminosity classification 



' For details see the Nightfall User Manual 

by Wichmann (1998) available at th e URL: 

|http://www.hs.uni-hamburg.de/DE/Ins/Per/Wichmann/Nightfall.htmI| 



Fig. 7. Variations of the mean, highest and lowest visible radius 
and temperature for model E5. Left: Values for the primary star 
Right: Values for the secondary star. 



problem. To refine our classification, we again used the Walborn 
& Fitzpatrick (1 19901 1 atlas. This atlas was used with the indi- 
vidual spectra, which for real observational data can only be as- 
sessed through spectral disentangling. The atlas is more qualita- 
tive than the Conti-Mathys criterion but relies on a large number 
of lines. We can now achieve a better agreement with the classi- 
fication of the real binary systems (see Table|5). 

4.2. Struve-Sahade effect 

In our sample of binaries, two are known to present this effect: 
HD 93403 and HD 152248. They inspired our models E2 and 
E5. In our analysis, we detected variations of the line strength in 
models El, E2, E4, and E5. However, we observe an S-S effect 
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Table 4. Para meters ( from Antokhina et al. (EOOOt . Rauw et al. (ISOOOt . Gies et al. (|2002] |. Sana et al. (l2008l l. Mayer et al. (l2008l l. 
and Sana et al. j2001l) for models El to E5) of the eccentric binary systems. Inclinations in brackets stand for non-eclipsing systems. 



Parameters 


Model El 


Model E2 


Model E3 


Model E4 


Model E5 


Period (day) 


6.08 


15.093 


9.6466 


5.604 


5.816 


Eccentricity 


0.46 


0.234 


0.156 


0.259 


0.133 


Argument of periastron of secondary (°) 


197.4 


202.5 


302 


284 


82 


Mass ratio 


0.44 


0.54 


0.63 


0.74 


0.98 


Inclination (°) 


(60) 


(32) 


(56) 


(60) 


67.2 


Mass of primary (Mq) 


45 


68.5 


36.2 


24.5 


27.8 


Mass of secondary (Mq) 


20 


37.3 


22.9 


18.2 


27.2 


Primary polar temperature (K) 


49000 


39300 


40500 


32800 


34350 


Secondary polar temperature (K) 


36500 


40100 


35000 


31200 


34000 


Polar radius of primary (Rq) 


9.2 


24.0 


9.5 


10.3 


16.0 


Polar radius of secondary (Rq) 


7.0 


10.0 


7.7 


7.9 


14.5 


visin(i) (km s"') 


135 


144 


102 


152 


135 


V2Sin(i) (km s"') 


145 


75 


164 


133 


135 


jS*"' of primary 


0.64 


2.06 


1.80 


1.09 


0.81 


jS*"' of secondary 


0.91 


2.59 


3.57 


1.24 


0.89 


TIDES code parameters 


Viscosity, v, of primary (R^day^') 


0.05 


0.045 


0.028 


0.05 


0.05 


Viscosity, v, of secondary (iJgday"') 


0.05 


0.01 


0.028 


0.02 


0.02 


Layer depth 


0.01 


0.1 


0.07 


0.02 


0.1 


Polytropic index of primary 


1.5 


3 


1.5 


3 


3 


Polytropic index of secondary 


1.5 


1.5 


1.5 


1.5 


3 


Number of azimuthal x latitudinal partitions 


500 X 20 


500 X 20 


500 X 20 


500 X 20 


500 X 20 



Notes. '"'The /3 parameter measures the asynchronicity of the star at periastron and is defined by yS = 0.02^^ x j^-^^, where v,,,, is the equatorial 
rotation velocity, R is the equilibrium radius, and e is the eccentricity. 

Table 5. Spectral classification using the Conti-Mathys criterion and the Walborn & Fitzpatrick atlas. 



Stars 


Combined spectra 


Individual spectra 


W&F atlas 


Observational analysis 


Model El' 


03 


03 


03-4 (IV)-V 


03V 


Model EP 


07.51 


07I 


08-8.5 III-(V) 


08V 


Model E2' 


05.5 


04 


05-6 (I)-III 


05.5I 


Model E2^ 


04 


04 


07-8 III-(V) 


07V 


Model E3' 


061 


05.5I 


06.5-7 III-(V) 


06.5V 


Model E32 


06.51 


07I 


08.5 V 


08.5V 


Model E4' 


08I 


08.5III 


09-9.5 (III)-V 


09III 


Model E42 


09III 


09.5V 


09.5-9.7 III-(V) 


09.7V 


Model E5' 


06.5I 


07I 


07-8 III-(V) 


07III 


Model E5- 


07I 


07I 


07.5-8.5 (III)-V 


07.5III 



Notes. 'Primary. ^Secondary. 



in the sense of Linder et al. ( I2007I I only for some lines in model 
E2. 

For all models, we investigated eight "well-chosen" lines in 
the spectra at 20 phases (the phase zero corresponds to the peri- 
astron passage). The lines were chosen following three criteria: 
first, we selected the lines for which an S-S effect is reported 
in the literature. The second criterion is based on a visual de- 
tection of variations in the line profile in the synthetic spectra, 
and finally, the third criterion is based on the fact that some lines 
present the S-S effect more often than others. As in Paper I, we 
draw attention to possible blends with nearby lines that can mod- 
ify the strength and lead to misinterpretations. 

We followed the same procedure as in Paper I and measured 
the EWs in the simulated combined spectra of the binary at dif- 
ferent phases as well as in the simulated spectra of the individual 
components of the binary at the same phases. This approach al- 
lows us to compare the line strengths that are deblended from 
the combined spectrum with the actual individual spectra, which 
is not possible with the real observational data. We found varia- 



tions in the EWs in many of the individual spectra as well as in 
the combined spectra. However, the agreement between the re- 
sults of the deblend/line routine and the actual individual spectra 
is mostly poor We present in Figs. |9] and [TOl examples of typi- 
cal variations observed in our analysis. The (inverse) U-pattern 
(Fig.^right, grey or FiglTO]Ze/f, black) is often observed and is 
directly related to the orbital motion. Constant EWs or irregular 
variations are also often observed. A constant EW is mostly ob- 
served in long-period systems. We also encountered other types 
of variation patterns such as those presented in FigJTO] (Wg/ir, 
grey) during our analysis. The amplitude of the EW variations 
depends on the line and the system. The variations observed on 
individual spectra can be explained by the physical variations 
of the temperature and gravity at the stellar surface during the 
orbital cycle (we give an example of the interpretation of the 
variation in subsection l4.2.1b . The explanation of the S-S effect 
that we suggested in Paper I seems to be reinforced by these 
new analyses. The S-S effect seems to be due to the combina- 
tion of the spectra in which the lines have a non-Gaussian/non- 
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He I A 4471 



Hell A 5412 



- Primary 

- Secondary 

- Combined 

- Midas fit 




■ -* \ 
1 \ 



5580 55S5 5590 5595 5600 5605 

Wavelength (A) 

Fig. 8. Combination of the O III /15592 line of the primary and 
the secondary for model E2 at phase 0.25. The deblend/line rou- 
tine will not fit the line of the secondary properly. Top: The re- 
sulting lines of the two Gaussian fits. Middle: Individual lines 
of the primary and secondary computed with our model and nor- 
malized to the continuum of the entire system. Bottom: Spectrum 
of the system and corresponding fit of MIDAS. 



Lorentzian and even asymmetric profile. Under these conditions, 
the deblend/line routine introduces systematic errors that mimic 
the S-S effect. This is illustrated in Fig.|8] which shows that the 
deblending routine erroneously selects the bump that is created 
by the superposition of the two lines as the secondary star's ab- 
sorption line. The secondary's line appears as a weak contribu- 
tion on the red wing of the primary star's absorption. However, a 
fit with two Gaussian profiles will consider the bump as the core 
of the secondary line and thereby overestimate its strength. 



He I A 4471 



O III A 5592 




0.95 f % 
0.85 « • 





0.5 1 

Orbital phases 
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Fig. 9. Example of EW variations for the He I /I 447 1 line for 
model E2 and for the OIII A 5592 line for model E4. The 
EWs have been normalized to the highest value. Black: Primary. 
Grey: Secondary. Top: Variations measured on individual spec- 
tra (EWp, max, Hel = 0. 1 7 A, EWp, max, OIII = 0. 1 8 A, EW^, Hel = 

0.34A, EWs, max, OIII = 0.11 A). Bottom: Variations measured on 
combined spectra (EWp_ max, Hel = 0.17A, EWp,max,oiii = 0.14A, 

EWs, max, Hel = O.lOA, EWs. max, OIII = 0.08A). 




0.5 

Orbital phases 
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Fig. 10. Example of EW variations for the He I /I 447 1 and 
Hell A 5412 lines for model El. The EWs have been nor- 
malized to the highest value. Black: Primary. Grey: Secondary. 
Top: Variations measured on individual spectra (EWp. max, Hei = 

0.04A, EWp,max.HeII^ = 1.07 A, EWs, max, Hel = 0.59A, 

EWs, max. Hell - 0.98A). Bottom: Variations measured on com- 



bined spectra (EWp,max,HeI 
EWs, max, Hel - 0.21 A, EWj, max. Hell 



0.05A, EWp 

, max, Hell 

0.29A). 



0.80A, 



4.2.1. Model E1 

In this model, we mainly observed variations illustrated in 
FigHni These variations can be explained by the orbital motion 
of the stars and the variation of the visible part of the stars. At pe- 
riastron passage, we see the rear side of the secondary. Because 
the stars are close, they are quite deformed and therefore the 
visible part of the secondary is cool, which leads to a reinforce- 
ment of the He I lines. For the primary the hemisphere facing 
the companion is visible at this phase. However, because the pri- 
mary is deformed, and owing to the large temperature difference 
between the two stars, the reflection process does not heat the 
primary significantly. Thus it is also the coldest part of the pri- 
mary that is visible at this phase. When the binary separation 
increases, the stars become nearly spherical. The temperature of 
the primary becomes nearly constant and the EW of its lines 
do not change. However, for the secondary, we begin to see the 
front part of the star Whilst the reflection is not very eff'ective 
for heating the primary, it is very effective for the secondary and 
accordingly, when we see the front part of the latter, the temper- 
ature increases. The increase is stronger because the separation 
between the stars decreases during the second half of the orbital 
cycle. 

The EWs measured on combined spectra (see FigJTOl 
Bottom) are not consistent with those measured on the individ- 
ual spectra (except for the Hel /i4143 line of the secondary). 
As mentioned above, this inconsistency is due to the deblend- 
ing routine, which does not properly separate the lines of the 
primary and secondary stars. For example, the blue-shifted com- 
ponent may have a larger EW than when it is red-shifted. This 
inconsistency is reduced when the lines are measured at orbital 
phases at which they are clearly separated. 
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Fig. 11. Example of EW variations for He I A 5016 and Hell 
A 4200 for model E5. The EWs have been normalized to the 
highest value. Top: Variations measured on individual spectra 

(EWp,max,HeI = 0.20A, EWp, ^a^. Hell = 0.48A, EWj, max, Hel = 

0.21 A, EWs. max, Hell = 0.48A. Bottonv. Variations measured on 
combined spectra (EWp_ max, Hei - 0. 12A, EWp_ max, Heii = 0.27 A, 

EWs,max,HeI = 0.11 A, EW^, max. Hell = 0.24A). 



4.2.2. Model E5 

For this last model, we investigated the He I AAAQlb, 4143, 
4471, 4713, 5016, SilV ^4089, Hell .14200, 4542, and OlII 
A 5592 lines. The EWs of the He I and Si IV lines in the primary 
display a U-pattern with a maximum at phase 0.0 - 0.05 and a 
minimum at phase 0.45 - 0.55. The EWs of the secondary dis- 
play a symmetric variation with respect to the primary, however, 
the amplitude of the variation is smaller The O III line displays 
the inverse behaviour compared to the previous lines. Finally, 
the He II lines in the primary display two maxima at phases 0.3 
and 0.7 with an intermediate minimum at phase 0.5 that give us 
an M-pattem variation. The secondary presents a low-amplitude 
U-pattern variation. 

The EWs measured on the combined spectra are nearly con- 
stant for the He I AA 4471, 5016, Hell /1 4200, 4542, and 
O III A 5592 lines for the primary and secondary stars. The Si IV 
line displays irregular variations but a mean value of the primary 
higher than the secondary for phases before 0.5 and a mean value 
of the secondary higher for the second part of the orbital cycle. 
The He I /14713 line displays constant EWs for phases before 
and after 0.5 but the value is different before and after phase 
0.5. Before phase 0.5 the primary has a lower EW than the sec- 
ondary and after phase 0.5 the primary has a larger EW than the 
secondary. Finally, the variations observed for the He I /1 4143 
line are consistent with the variations measured on individual 
spectra. All variations measured in our model underestimate the 
observed S-S effect. 

In conclusion, our models display, in general, some varia- 
tions of EWs during the orbital cycle. The variations on the spec- 
tra of individual components are caused by the orbital motion 
and the modulation of the visible part of the star. However, the 
variations that we measure on combined spectra and the visual 
variation of the relative strength of the line of the primary and 
secondary stars are not strong enough to completely explain the 
S-S effect as observed by Sana et al. (i200ri l. 
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Fig. 12. Synthetic light curves of model E5 in the range AA 
3500, 7100. The light curves have been computed with the 
CoMBiSpeC + TIDES models (black) and with the Nightfall 
model {grey). 



4.3. Synthetic broad-band light curves and radial velocity 
curve 

The synthetic light curves were computed in the wavelength 
range 3800 - 7100 A. We compared our light curves to those 
obtained with Nightfall. Nightfall provides synthetic light curves 
based on an instantaneous Roche potential and thus does not take 
any viscous stress into account. We point out here the similar- 
ity and difference of the two models. For models El and E4, 
our light curves present qualitatively the same characteristics as 
Nightfall. The TIDES code leads to more compact objects than 
the instantaneous Roche potential. Therefore, because the stars 
are more compact, all other parameters being equal, the gravity 
darkening is less strong and the reflection effect is more effective 
and consequently the stars are hotter Near periastron, the differ- 
ence between the two models is stronger because our model pre- 
dicts a weak variation of 0.02 mag and Nightfall predicts no vari- 
ation. The very low inclination (60°) explains this small change 
of magnitude. Because the stars are less deformed (than in the 
Nightfall model), the reflection effect is more effective near pe- 
riastron and accordingly, it also increases the temperature of the 
stars which enhances the difference between the two models. 
Some other differences arise because of the time delay induced 
by viscous stress. For model El, the agreement between the two 
models and the light curve observed for HD 93205 by Antokhina 
et al. (12000) is poor, however Models E2 and E3 give very dif- 
ferent results with our code and Nightfall. These two systems 
have a longer period and are more compact in our model than 
in Nightfall. The discrepancies probably come from the different 
radii of the stars. Finally, the light curve of model E5 displays 
the two eclipses in both models but in our model, the eclipses 
are deeper than evaluated with Nightfall (see FigfT2ll. Again the 
difference can be explained by the smaller radius of the star in 
our model. Some other small differences can be easily explained 
by the viscous effects. The agreement between our model and 
the observations of Mayer et al. (2008 l is rather good, though 
we slightly overestimate the depth of the eclipse of the primary. 

The mean brightness ratios in the wavelength range 3800 - 
7100 A are similar to the literature values. Finally, we deter- 
mined the amplitude of the radial velocity curves and found a 
good agreement with the observations (see Table |6) for all sim- 
ulated systems. 



5. Summary and perspectives 

We have presented improvements of our mathematical model 
that allow us to compute the physical properties on the surface 
of massive stars in binary systems that contain main-sequence 
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Table 6. Semi-amplitude of radial velocity curves for eccentric 
models and comparison with observations (in kms"') 



Stars 


K, 


K, 


Kl_obs 


Kl, obs 


Model El 


136.3 


308.8 


132.6 


313.6 


Model E2 


77.3 


142.8 


79.3 


139.0 


Model E3 


125.3 


197.5 


117.0 


211.0 


Model E4 


157.6 


213.9 


162.4 


213.9 


Model E5 


209.8 


208.4 


206.9 


211.7 



0(B) stars. The first improvement is the inclusion of the radi- 
ation pressure effect on the shape of the stars. The second im- 
provement is the use of the TIDES code to compute the shape 
of the stars in eccentric and/or asynchronous systems. In both 
cases, we took into account various effects such as gravity dark- 
ening, reflection, and limb-darkening which allowed us to com- 
pute the temperature distribution at the stellar surface. Then we 
used the TLUSTY OSTAR2002 and BSTAR2006 grids (Lanz & 
Hubeny !2003l 1200 7 ) to compute the spectra of each star of the 
system as a function of orbital phase. We showed that in a num- 
ber of cases, the radiation pressure does not have a strong impact 
on the shape of the stars and therefore on the spectra for the mod- 
els studied in Paper I. This implies that the conclusions of Paper I 
for the S-S effect remain valid. Our results showed that the radia- 
tion pressure has a weak impact on the shape of highly deformed 
stars of (over)contact binary systems like model 3. However, in 
model 3, the small changes in the shape of the stars have resolved 
our previous problem of the surface temperature distribution that 
we failed to reproduce with our first version of the algorithm . We 
also studied the impact of the gravity-darkening parameter vari- 
ation on the spectra of model 3. Our results indicate that, given 
the assumptions made for the computation, there is no need to 
use a more complex treatment of the gravity darkening than the 
classical von Zeipel theorem. The model 4 clearly showed that 
radiation pressure can have a strong impact on the shape of the 
stars. Many of the spectral lines in this model display phase- 
locked profile and/or strength variations. However, these varia- 
tions are not sufficient to reproduce the S-S effect observed by 
Bonanos ( |2009l ). 

The second part of this paper investigated the eccentric sys- 
tems. We studied many lines and detected phase-locked profiles 
and/or strength variations in many cases in the individual spectra 
of both components. These variations are caused by the change 
of orientations of the stars as a function of orbital phase. They 
thus reflect the non-uniform temperature distribution across the 
stellar surface. The variations measured on the combined spec- 
tra of the systems, however, often disagree with the measure- 
ments made on individual spectra. Our results also showed that 
the deblending routines that fit two Gaussian or Lorentzian pro- 
files often fail to measure the line EWs properly and lead to in- 
correct interpretations. The variations measured in models E2 
and E5 underestimate the S-S effect observed in HD 93403 and 
HD 152248, which in turn might be caused by our not predict 
intrinsic asymmetric line profiles in these cases. In Paper I, this 
asymmetry was our explanation for the S-S effect: the sum of 
asymmetric intrinsic line profiles of individual stars are not prop- 
erly deblended by Gaussian/Lorentzian profile fits. The absence 
of this asymmetry in eccentric systems is, however, not easy to 
explain and could be linked to weaker interactions between the 
stars in eccentric binaries than in close circular ones. Therefore, 
we cannot generalize the interpretation of the S-S effect of Paper 
I to the eccentric systems studied here. However, this highlights 
that analyses of binary system spectra are very difficult and some 



observational effects could be generated by analysing techniques 
and not by physical processes. 

For clearly detached systems, spectral disentangling 
(Hadrava [T9951 Gonzalez & Levato l2U06l Simon & Sturm 
119941) helps to overcome the difficulties of the deblending 
routines. However, for systems with a strong temperature 
gradient at their surface (e.g. HD 100213, Linder et al. l2007l l. 
the technique fails because different lines have different radial 
velocity amplitudes. In eclipsing binaries, disentangling cannot 
be used for the phases near the eclipses. This technique provides 
mean spectra of both components, and even if we can measure 
variations with respect to this mean, we cannot compute the 
spectra at each orbital phase. 

A future step could be to introduce cross-talk which is surely 
substantial in eccentric and asynchronous systems. Finally, we 
plan to include a wind interaction zone between the stars which 
could contribute to the heating of the stellar surface in two ways, 
either by backscattering of the photospheric photons, or by ir- 
radiation of X-ray photons emitted by the shock-heated plasma 
in the wind interaction zone. This requires, however, more so- 
phisticated atmosphere models than what we have employed so 
far 
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